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The coupling qe'e-k is related to the form factor at zero momentum of the axial current between B* and B 
states. Moreover it is related to the effective coupling between heavy mesons and pions that appear the heavy 
meson chiral Lagrangian. This coupling has been evaluated on the lattice using static heavy quarks and light quark 
propagators determined by a stochastic inversion of the fermionic bilinear. We found the value g = 0.42(4)(8). 
Beside its theoretical interest, this quantity has phenomenological implications in B — > tt + decays. 



1. INTRODUCTION 

A matrix element with important phenomeno- 
logical applications is 

{B\p)7:+{q)\B*+{p')) (1) 

which can be parametrised in terms of the form 
factor ggg,^{q^) in the following way: 

- 5B*s.(g')g^?7^(p')(2^)''5'(p' -p'-q) (2) 

where rj^ij)') is the polarization vector of the 
asymptotic state B*^{p'). 

For on-shell external states, p' = p + q, one can 
perform an LSZ reduction of the pion in eq. 

i{ml ~ q') I e»«- {B{p)Mx) \B*{p')) d^x (3) 



Using the PCAC definition of the pion in terms 
of the axial current 

one obtains 



(4) 



, e^^-{B%p)\A^{x)\B*{p'))6?x (5) 

Since in the limit g ^ 0, eq. (||) and eq. (0) 
parametrise the same matrix element of eq. (O, 
one derives a Goldberger-Treiman like relation for 
the pion system 
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where g is defined as 
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The parameter g is also the effective coupling 
which appears in the the Heavy Meson Chiral 
Langrangian Q and can be used to constrain the 
form factors that appear in the matrix elements 
of the weak vector current, u'y'^b. Which matrix 
elements are relevant B ^ tt + II decays. One 
of the best present techniques to extract Kib is by 
comparing the experimental values for these form 
factors with the theoretical ones, obtained by fit- 
ting lattice Montecarlo results. Therefore a pre- 
cise determination of g can be used to reduce the 
number of free parameters in the fit and, there- 
fore, reduce the theoretical uncertainty on Uub- 

Even if the matrix element of eq. (|]) does not 
directly correspond to a physical process, because 
B* —>■ B+TT is kinematically forbidden, the equiv- 
alent process for D systems, D* D + tt, is in- 
deed physical and the effective coupling for this 
D decay process coincides with our g up to 1 /ttIc 
corrections. 

2. SIMULATION 

We have performed our simulations on 20 
quenched gauge configurations, generated on a 
12^ X 24 lattice at /3 = 5.7, corresponding to 
a^^ = 1.10 GeV, with a tadpole improved SW 
action {csw — 1-57). The heavy quark propaga- 
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tors are evaluated in the static limit, while the 
light quark propagators are evaluated performing 
a stochastic inversion on the fermion matrix Q 

iQ-')^j = J [d0](Q,fe0fc)*0,e-^'(«''3).™0™ (8) 

and 10 pseudofermionic fields 0i for each gauge 
configuration were generated by Montecarlo. 
Moreover the maximal variance reduction tech- 
nique is implemented to reduce the statistical 
noise Two values of k are considered, ki — 
0.13843 and K2 = 0.14077, corresponding to a 
light quark mass of 140 MeV and 75 MeV respec- 
tively, with a critical value Kcrit — 0.14351. The 
interpolating operator jt (J) that create (annihi- 
late) a static B meson is smeared with a two-steps 
fuzzing procedure for the light quark fields. We 
indicate the local operators with the superscript 
L and the fuzzed ones with the superscript F. 

We have evaluated the two point correlation 
function at zero momentum 

^FFiFL) ^ (0|J(y,0)Jt(y,<)|0) (9) 

both for two- fuzzed (FF) and one-fuzzed one-local 
(FL) J operators. The average is on the spa- 
tial position y and it has the effect of increas- 
ing the statistics. This has been possible thanks 
to all-to-all stochastic propagators. From fitting 
the asymptotic (large t) behavour of C2 with a 
single exponential we have extracted Z^^^^ = 
(0| J^(^) \B) /V27^. 

Analogously we have computed the three point 
correlation function at zero momentum 

(0| J(y,-ii)A^(x + y,0)Jt(y,t2)|0) (10) 

for two-fuzzed J operators. Here the average is on 
both the spatial coordinates y and x, but keeping 
fixed r — |x| . For the spatial indices n = 1, 2, 3 we 
also used the rotational invariance of the lattice 
so that C31 = C32 = C33. 
Finally we have computed 

E^r,t) = iZ^rci/{r,t,t)/[cr{tf (11) 
and 

E{r,t) = {Z^fCif[r,t,t)/ [C|^^(i)]' (12) 



These two functions of r are plotted in fig. |l| for 
the lightest value of k and values oit = 3, 4, 5, 6. 
As expected E^{r,t) ~ 0, because the polariza- 
tion of a static B* meson is ortogonal to the tem- 
poral direction, while the averaged spatial com- 
ponent E[r, t) is almost independent on t and 
presents an exponential decay in r (this is con- 
firmed by a number of different fitting tests and 
by visualizing the points in log scale). In fact in 
the asymptotic regime in t 

E{r,t) = {B\A^[r)\B*) /(2mB) (13) 

Going back to the definition, eq. (0), the lattice 
regularised g coupling can be extracted from the 
spatial integral of /(r) = Se~'^/^°, the function 
fitting 'E(r,t). We find 

V... = //„(4„-.,d.- = {J;^^g (14) 

We observe that the main source of error in g'*^** 
is the lattice breaking of rotational invariance, in 
fact the distribution of the points E{r) around 
the fitting function /(r) exhibits a regular pattern 
which is independent on t and k. 

3. MATCHING 

From C2 we are able to extract the B meson de- 
cay constant (in the static approximation), using 
the relation 

^static ^ ^static /^^L ^-3/2 

V rriB 

For our lattice — 0.78 is computed using 

the Lepage-Mackenzie procedure of matching 
lattice vs continuum at a best scale, which in our 
case is determined to be q*a = 2.29. Our result is 
j^tatic ^ o.43(l)(8) GeV. We observe that a one- 
mass fit gives a value of f^'^*'^^ 20% bigger than 
the value extracted from a three-masses fit . In 
fact we obtain a value for this quantity that lies 
above other lattice calculations of the same quan- 
tity. We take into account this effect by adding 
a systematic error to our results (in particular to 
g) of the order 20%. 

The quantity we are interested in is the g cou- 
pling renormalized in the MS scheme at the ms 
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E"{r, t),K^ 0.14077 




4 6 
distance r 



E{r, t),K = 0.14077 




4 6 
distance r 



Performing a naive extrapolation to the chiral 
limit and including our evaluation for the system- 
atic error we summarise our result as 



g = 0.42(4)(8) 



(18) 



This number has to be compared to the best esti- 
mate from a global analysis of available results ||] 
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0.38 



(19) 



and with an estimate obtained from lattice data 
for the semileptonic B decay form factor (assum- 
ing Vector Meson Dominance), 



g = 0.50(5)(10) 



(20) 



A recent, independent, lattice analysis con- 
firmed our result 

One remark is in order. We consider this study 
an exploratory one because of the small lattice, 
the large lattice spacing, the poor chiral extrap- 
olation and the use of a non-standard pioneer- 
ing technique such as stochastic propagators. As 
we observed, the g parameter has important phe- 
nomenological applications and we believe our 
study shows that this number can be evaluted 
on the lattice. A better determination is required 
and is feasible with present computing power. 



Figure 1. Plots of i?°(r, t) and E{r, t) as funtions 
of r for different values of t {— 3, 4, 5, 6) for H2. 
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(16) 



^0 



where Z^'°°p = 1 - 13.8^, u},-'°°^ = 1 



TT^^f and uq — 0.86081 is the average plaquette. 
These values imply z^^'^^°^° = 0.806. 
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4. CONCLUSIONS 

Substituting the values of eq. 
we obtain 

_ f 0.44(4) for Ki 
^ ~ 1 0.43(4) for K2 



§ into eq. (16) 
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